Designer switches: Effect of crystal planes on time-dependent electron transport 

through an interacting quantum dot 
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The time-dependent non-crossing approximation is utilized to determine the effects of the crystal 
planes of gold contacts on time dependent current through a quantum dot suddenly shifted into the 
Kondo regime via a gate voltage. For an asymmetrically coupled system, instantaneous conductance 
exhibits complex fluctuations. We identify the frequencies participating in these fluctuations and 
they turn out to be proportional to the separation between the sharp features in the density of states 
and the Fermi level in agreement with previous studies. Based on this observation, we predict that 
using different crystal planes as electrodes would give rise to drastically different transient currents 
which can be accessed with ultrafast pump-probe techniques. 
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PACS numbers: 72.15.Qm, 85.35.-p, 71.15.Mb 

Time-dependent electron transport in single electron 
devices is a subject of fundamental importance in molec- 
ular electronics since it is widely believed that these de- 
vices have the potential to replace the conventional MOS- 
FET transistors 1 - in the near future thanks to the rapid 
advances in nanotechnology. Detection of electrons in 
real time 2 - is expected to play an important role in devel- 
opment of quantum computers^ and single electron guns 4 
as well. 

Time-dependent current arising from sudden switching 
of the gate or bias voltage^— has been shown to exhibit 
various time scales 8 -* 9 -. Moreover, interference between 
the Kondo resonance and the sharp features in the con- 
tacts' density of states emerges in the long timescale asso- 
ciated with the formation of this many-body resonance 1 ^. 
Effective one-electron theories indicate that the transport 
properties in steady state depend on the the type of elec- 
trode metalii, the contact-structure^ 2 - and the indices of 
crystal planes of electrode metali^. Same approach also 
predicts that the electrode metals would alter the recti- 
fying performance of the device^ 4 -. 

In previous studies of single electron devices, both the 
Green's function technique o 10 ' 15 ' 16 and nonequilibrium 
diagrammatic Monte Carlo metho d 17 ' 18 showed that the 
band structure of the contacts has significant influence 
on the shape of the transient current. In all these cases, 
simple unrealistic bands have been used for ease in cal- 
culations. In this letter, we will carry out a compar- 
ative study using three different crystal planes of gold 
as contacts to determine the transient current through a 
single electron device in the Kondo regime. A unique ap- 
proach which involves using the outcome of an ab-initio 
calculation as an input in a many-body technique will be 
adopted. 

Single impurity Anderson Hamiltonian representing a 
discrete spin degenerate level of energy Cdot connected to 
Fermi liquid electrodes describes the physical behaviour 



of this system adequately. Upon applying the auxiliary 
boson transformation to the Anderson Hamiltonian, it is 
transformed into 

H{t) = ^2\e k n kaa + V a (e ka ,t)cl arj tff a + R.c. 

kotcr 

y^e<fot(*)rc<r) (1) 
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where operators fl{f a ) and c\ arj {c kaa ) with a=L,R cre- 
ate (destroy) an electron with spin a within the dot and in 
the left(L) and right(R) contacts respectively. The corre- 
sponding number operators are n a and n kaa whereas the 
hopping amplitude for each contact is denoted with V a . 
Operator b'(b) creates (destroys) a massless boson within 
the dot and double occupancy is prevented by setting the 
sum of the number of bosons and the electrons to unity. 
If one neglects the explicit time dependency of the hop- 
ping matrix elements, the coupling of the quantum dot 
to the contacts can be cast as T L ^(e) — Tl(r)^l(r)(^)- 
In this expression, ^l(r) is a constant determined by 
f L(R) = 2n\V L ( R )(ef)\ 2 and £ L (fl)0) is the density of 
states of the contacts. 

The density-functional theory (DFT) calculations for 
the density of states (DOS) of the Au contacts have 
been performed with the full-potential linearized aug- 
mented plane wave (FP-LAPW) method using WIEN2K 
package 19 . As shown in Fig. [TJ three geometries, i.e. 
(OOl)-surface, (lll)-surface and (lll)-pyramid, have 
been adopted to simulate the Au electrode profile. In 
order to build the structures of these three geome- 
tries, atomic slabs with 13 Au layers have been used. 
The distance between the two opposite electrodes is 30 
Bohr radius in all three geometries. The exchange- 
correlation potential of the generalized gradient approxi- 
mation within the Perdew, Burke, and Ernzerhof (GGA- 
PBE) form- 2 - has been used in all calculations. The plane 



2 



(a) 



. * . . • 

to WW to« 



• V to fete «r* 

• to-* VW* 



(C) 



to V 



« to to 



» ». » i 
w to to 



t : 



to v 

to * to 

to to 



to to 



> ' 



FIG. 1: This figure shows the three geometries, (a) (001)- 
surface, (b) (lll)-surface and (c) (lll)-pyramide, used to sim- 
ulate the profile of Au electrode. The Au atoms in blue color 
are the atoms relevant to the transport properties between 
electrodes in all three structures. 



wave cut-off has been determined by R m tK max = 6.5 and 
Imax — 10. K-mesh of 30x30x3, 36x36x2 and 6x6x2 have 
been adopted for (OOl)-surface, (lll)-surface and (111)- 
pyramid respectively. Note that only the DOS of topmost 
atoms, which are relevant to the transport properties be- 
tween electrodes in the surface/pyramid structures, will 
be used in the following many body calculations. The 
resulting density of states is shown in Fig. [2j 

We performed a fitting procedure involving a linear 
combination of Gaussians given by 



p(«) = £ 



exp (-2(—^n (2) 



to the actual DFT data such that all of the sharp peaks 
would be captured. An extra broad Gaussian has been 
added in each case to ensure that the density of states at 
the Fermi level is non-zero and the entire bandwidth of 
the material is covered. The outcome of this data fitting 
procedure that involves Gaussians of varying linewidth 
and peak position is also shown in Fig. [2] for three differ- 
ent geometries. In the remaining part of this letter, we 
will shift to atomic units where h = kg = e = 1. 

The Dyson equations are solved for the retarded and 
less than Green's functions in a discrete two-dimensional 
grid. This requires the self-energy of electron and mass- 
less boson as an input. We use non-crossing approxi- 
mation to define the self-energies and close these cou- 
pled integro-differential equations^^. The net current 
follows from the resulting Green's functions denoted by 
Gpseu {t,t') and £?<( fl ) (t, f'). The general expression for 
the net current I(t)22- can be recast by using the slave bo- 
son decomposition method^. The equation we will use 
in this letter to determine the net current is then written 
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FIG. 2: Density of states of (001) surface, (111) surface 
and (111) pyramid calculated using DFT is shown with black 
dashed curve from top to bottom respectively as a function 
of separation from the Fermi level. Each Gaussian used to 
capture the sharp features and the Fermi level is also shown 
with green curves in each geometry. 



as 



J(t) = -2(f L ~ f R ) Re / dt^ {t, tx)h{t - h) 



2T L Re[ / dti(UMi)+£ u (Mi))/ L (t-ti) ~ 



2T R Re[ / dfci(&(-,ti) + £„(*, ti))/fl(*-*i) 



(3) 



where UMi) = G< seu {tM)B R {ti,t) and £u(Mi) 
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G* eu (t,t 1 )B<(t 1 ,t). InEq. ©, f L (t-tx) and J R (t-ti) 
are the convolution of £(e) with the Fermi-Dirac distri- 
butions of contacts and h(t — ti) is the Fourier transform 
of £ (e)^£. The conductance G is equal to the current di- 
vided by the bias voltage V . In the following discussion 
f] = jr^, where Ft of = Tl + f r, will be referred as the 
asymetry factor. 

Kondo effect is a many-body resonance pinned to the 
Fermi levels of the contacts in the dot density of states 
and manifests itself as an enhancement in the conduc- 
tance at low temperatures. The linewidth of the Kondo 
resonance is characterized by an energy scale Tk (Kondo 
temperature) given by 

(DT tot \^ ( ^|e dot |\ 

Tk * {—) exp lr V) ■ ( 4 ) 

where D is the half bandwidth of the conduction electrons 
and Ttot = f£(e/)- 

We will investigate the transient current for the case 
where the dot level is switched from t\ = —4T to t to 
£2 = — 2r tot at t=0 with a gate voltage. For all three ge- 
ometries, a transition from a non-Kondo state to a Kondo 
state takes place. Note that the Kondo temperature in 
the final state would be slightly lower for (111) pyramide 
than the other two geometries as a result of shorter con- 
duction electron bandwidth as seen in Fig. [2] 

Transient current in Kondo timescale is depicted for 
all geometries in Fig. |3] after switching to the final dot 
level in infinitesimal bias. We do not display the short 
timescale where where the conductance reaches a maxi- 
mum before it starts to fall off for large asymmetry fac- 
tors due to striking similarity with previous studies^. 
The most remarkable feature is the drastic difference in 
conductance fluctuations. This effect stems purely from 
the difference in band structure of the contacts since all 
other parameters are held constant. We should point out 
that the slight difference between the Kondo temperature 
of (111) pyramide and the other two geometries does not 
change the relative behaviour of the fluctuations for these 
geometries because T/Tk is known to alter only the am- 
plitude of the fluctuations 10 , not the overall pattern of 
the transient current. 

It is clear that these fluctuations are the result of an ad- 
mixture of sinusoidal oscillations with different frequen- 
cies and amplitudes. Each individual frequency can be 
extracted by taking the Fourier transform of the cur- 
rent. We found that the frequencies are proportional to 
the separation between the peak positions and the Fermi 
level for all cases. This explains why transient current for 
(001) surface exhibits a more erratic pattern compared 
to the others as five distinct frequencies are involved in 
it. Moreover, there are other peaks that appear in actual 
DFT data but we did not include them in our fitting 
since they have negligible contribution to the fluctuation 
pattern. The oscillation amplitude associated with them 
is probably too small to have a discernible effect because 
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FIG. 3: Panels from top to bottom show the instantaneous 
conductance versus time in Kondo timescale after the dot level 
has been switched to its final position for (001) surface, (111) 
surface and (111) pyramid respectively with an asymmetry 
factor of 0.9 at T=0.009rt ot in infinitesimal bias. 



the peaks are not prominent compared to the surround- 
ing structure. 

The results presented here describes a way to probe the 
detailed band structure of a molecular switch by measur- 
ing the transient current flowing through it. It is possible 
to capture the transient current fluctuations occurring 
in the femtosecond timescale with ultrafast pump-probe 
techniques^. This would pave the way for designing cus- 
tom switches for future organic computers. 
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